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Abstract 

I present an analysis of fully differential single-top-quark production plus jets at next-to-leading 
order. I describe the effects of jet definitions, top-quark mass, and higher orders on the shapes 
and normalizations of the kinematic distributions, and quantify all theoretical uncertainties. I 
explain how to interpret next-to-leading-order jet calculations, and compare them to showering 
event generators. Using the program ZTOP, I show that HERWIG and PYTHIA significantly 
underestimate both s-channel and t-channel single-top-quark production, and propose a scheme to 
match the relevant samples to the next-to- leading-order predictions. 
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I. INTRODUCTION 



The discovery of single-top-quark production will be the flagship measurement of run II 
at the Fermilab Tevatron. Preliminary run II results ^| are already surpassing the limits 
set by the CDF 3| and DO Collaborations using data from run I. Over the next few 
years, the single-top-quark cross sections will be accurately measured at both the Fermilab 
Tevatron 4] and the CERN Large Hadron Collider (LHC) [5]. In order to compare the 
coming measurements of single-top-quark production to theory, a detailed understanding of 
the predictions for jet distributions, and the associated uncertainties is required. 

The measurement of single-top-quark production will provide an excellent opportunity to 
study the charged-current weak-interaction of the top quark. Within the standard model, 
a measurement of the cross section allows direct extraction of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix element |V^t6p. In Fig.^we see that V^h appears in the leading-order 
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(LO) Feynman diag rams for t-channel production 
duction 3 ll^ ll^ ll^, and VTt-asso dated production 
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.^21 1^. s-channel pro- 
H iJH 21, Q- This 



paper focuses on s-channel and t-channel production, which have sizable cross sections, and 
may be distinguished experimentally by the number of b tags in the final state [2^. Since 
the lyt-asso dated production cross section is small at the Tevatron, and no differential 
next-to-leading-order calculation exists, Wt production is not addressed here. 



(a) 



(b) 



(c) 



FIG. 1: Representative leading-order Feynman diagrams for (a) t-channel, (b) s-channel, and (c) 
VFt-associated production of a single top quark. The CKM matrix element Vtb appears directly in 
the production diagrams. 



Models that extend the standard model often predict large corrections to s-channel or 
t-channel production, but not to both [24]. Anything that produces an anomalous coupling 
Q, H H 0, E , Q H 3 ) or flavor changing neutral current Q, Q , between the top 
quark and any other quark opens new production modes in the t-channel, but only slightly 
reduces the fraction of b jets in s-channel production. Conversely, any process that allows 
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a new intermediate particle, such as strong dynamics 
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, a charged top-pion 
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431 14|, would 



Kaluza-Klein modes of the W |41], or a W boson 
enhance s-channel production, but is highly mass suppressed in t-channel production. In 
fact, s-channel single-top-quark production has been found to be the most sensitive probe 
of high- mass charged vector W bosons 3| at the Tevatron [46] and LHC 47 1. 

Both s-channel and t-channel production are important in direct searches for CP violation 
-R-parity-conserving supersymmetry j5l|, and i?-parity- violating supersymmetry 
Both channels also contribute a significant background to all signals 




that include W + iets with, or without, b tags. These backgrounds appear in a number 

n n 

of Higgs search channels j57|] and other new physics, such as i?-parity-conserving j58|] and 
i?-parity-violating 0| supersymmetry searches. 

The analytic form of the fully differential next-to-leading-order (NLO) cross section for 
s-channel and t-channel production of a single top-quark plus one jet appears in Ref. 
The spin-averaged version of that calculation is available from this author in the form of a 
FORTRAN program called ZTOP, and is used in this paper. A recurring theme throughout 
the analysis of the kinematic distributions is that the next-to-leading-order calculation de- 
scribes jets, and not partons. This affects the questions that may be asked about the final 
state, and contributes additional challenges to the matching between theory and experiment. 

The ultimate goal of performing differential calculations at next-to-leading order is to 
improve the quantitative connection between theory and experiment. As important as the 
distributions themselves, is a reliable estimate of the uncertainties in these distributions. 
Previous estimates of the uncertainties in the parton luminosities and kinematics [2^ have 
been somewhat rough, and are revisited with more careful analysis. Analytic calculations 
are not typically compared to raw data, but rather to event generators, such as HERWIG 



and PYTHIA [61[ , that have been tuned to fit the data. It is vital to determine whether 
these event generators at least predict the shapes of the measured distributions accurately. 
We will see that the current event generators provide very poor predictions of both s-channel 
and t-channel production. 

This paper is organized as follows. A general search strategy for discovering single-top- 
quark production appears in Sec. II Al This lays the groundwork for discussing the event 
rates in terms of the number of jets in the final state. The next-to-leading-order distributions 
for both s-channel and t-channel single-top-quark production at the Tevatron are presented 
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in Sec. [H] The hazards of treating NLO calculations as parton calculations are explored, 
and subtleties in the interpretation of jet calculations are discussed in detail. The bottom- 
quark mass is shown to have no measurable effect on the distributions. All uncertainties in 
the inclusive and exclusive theoretical calculations are described in Sec. 11111 In Sec. IIVI the 



PYTHIA j61[ and HERWIG |60| event generators are shown to have poor representations of 
the measurable kinematic distributions. In Sec. IIV CI a simple scheme is proposed to match 
the event generators to the next-to-leading-order calculation. The paper concludes in Sec. |3 
with a few directions for improvement. 



A. Strategy for discovery 

The final state for both s- and t-channel production of single top quarks contains two jets 
plus the decay of aW boson. In order to extract the signal from the backgrounds, the decay 
of the W to an electron or muon, and a neutrino is used. The most serious background 
for single-top-quark production is the production of ti pairs j^]. The number of ti pairs 
that will be reconstructed in detectors at the Tevatron is about 5 times the total number of 
single-top-quark events in the 1 6-tag sample, and about 10 times the number in the 2 6-tag 
sample. Because there is already a top quark in the final state, the essence of any study will 
be to reduce the ti background, while maintaining a measurable signal. 

The most effective method to discover single-top-quark production uses a "jet veto," 
which was first proposed in Ref. ||23j. This method relies on the observation that in ti 
production, the extra top quark tends to produce additional high-i^^- objects (leptons or 
jets) in the final state. The basic strategy is to form a sample of two jets, a lepton, and 
missing transverse energy Ej-. If there are any additional isolated "hard" jets or leptons 
then the event is rejected. The term "hard" can be loosely defined to refer to any jet or 
lepton with Et > 15-25 GeV, and \ri\ < 2.5-4. Experimentally, the "jet veto" is identical 
to separating the final states into a ly + 2 jet exclusive sample and a + 3 jet inclusive 
sample. 

We observe in Fig. 2 of Ref. 0] that placing cuts anywhere in the region of E^j = 15-25 
GeV for every jet yields the same significance. The reason is that acceptance for events that 
contain an extra radiated jet in the signal increases as the cut is lowered (thereby failing to 
pass the cuts) at the same rate that the acceptance for the highest-ii^T jet increases. Hence, 



4 



in order to avoid a sharp drop in efficiency, the jet veto should be set at the same Et as the 
cut on the second-highest-i?T jet. 

The largest remaining background is Wjj production — dominantly Wcj + Wjj in the 
1 6-tag sample, and Wbb in the 2 6-tag sample. This background decreases as a function of 
reconstructed top-quark mass. Hence, a top-quark mass cut reduces this to an acceptable 
level. If a top-quark mass cut cannot be used, e.g., because of poor jet-energy resolution, 
other discriminates will have to be found. However, these will be ineffective without the 
corrections described in Sec. IIVI 

The detect,™ strategies of Refs. Q Q should be .e.,ed into a progra-n of f^ee sue- 
cessive measurements. The ffist stage measures the 14^ + 2 jet signal with at least 1 b tag. 
This sums over all single-top-quark production modes, and requires the least amount of 
integrated luminosity for discovery. As data is accumulated, the samples can be split into 
independent subsets with 1 or 2 6 tags. The 1 6-tag sample can be used to study t-channel 
production, and the 2 6-tag sample can be used to study s-channel production. 



II. NEXT-TO-LEADING-ORDER DISTRIBUTIONS 

Fully differential next-to-leading-order cross sections predict jet distributions; they do 
not predict parton distributions. The consequences of this can be subtle, and require care 
when comparing fixed-order calculations with data, or with showering event generators. To 
understand these comparisons, only distributions for top-quark production at the Tevatron 
are examined. Distributions for antitop-quark production at the Tevatron differ solely by 
the signs of the pseudorapidities, and figures would be redundant. Since the physics issues 
are identical, no jet distributions for the LHC are shown. 

The challenge of calculating differential cross sections is in how to deal with final-state 
infrared singularities. Different methods have different solutions in detail, but they all 
impose one of two conditions: include a jet definition, or include a hadronization function. 
The first consequence of including a jet definition in the base theoretical calculation is that 
the same choice of jet definition must be used in any comparison to data. This may appear 
obvious, but it contains a subtlety: how do you compare a jet that consists of a few partons 
with one that contains multiple particles? In Sec. Ill (^1 we will see the effect of the choice 
of jet definition on theoretical distributions. Comparisons between theory and experiment 



5 



are actually made at the event generator level, rather than directly with data. Hence, in 
Sec. IIVI we will adopt an effective solution for matching the NLO calculations to the event 
generators. 

Before presenting NLO distributions, we must first understand which distributions are 
appropriate to calculate. We might be tempted to look at the LO diagrams and ask for the 
corresponding NLO distributions, but this would not be useful in general. Consider the LO 
diagram for t-channel production in Fig. ^a). If we plot the NLO transverse momentum of 
the jet that includes the final-state d quark, we see in Fig. |21 that the calculation is not very 
stable. The massive dipole formalism (MDF) and phase space slicing (PSS) calculations of 
Ref. provide the same results at the same scale. The methods only differ in the bin that 
exactly includes pTd = GeV, which is an artifact of the different subtractions methods, 
and in any case not observable. 




10 20 30 40 50 60 70 80 
PTd (GeV) 



FIG. 2: Transverse momentum pxd of the "d"-jet in t-channel production at leading and next-to- 
leading order for the phase space slicing (PSS) and massive dipole formalism (MDF) calculations, 
and two scales, nit and the double deep- inelastic-scattering scales (Q^ and Q"^ + m1). Leading order 
includes a /C- factor of 1.09. 

The problem with Fig. |21 is that the shapes of the distributions are not stable compared 
to scale variation or perturbative order. At leading order the d quark recoils against the 
massive top quark, and so has a finite distribution at small pTd- At NLO the distribution 
takes an arbitrary value that depends on the choice of scale. Also the shapes at larger (i.e. 
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measurable) pxd change between orders, which leads to the false impression that the NLO 
jets are softer than the LO parton. This figure demonstrates the danger of interpreting jet 
results in partonic language. 

The distribution that should have been plotted is the transverse momentum of the highest- 
Pt jet. In Fig. |21 we see that this distribution is stable in shape at different orders in 
perturbation theory, and for different scales. In fact the difference between the LO and 
NLO distribution is smaller than the NLO scale uncertainty. We also see that the choice of 
scale effects only normalization. Hence, we have a very well-predicted quantity. 
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FIG. 3: Transverse momentum pxji of the highest-pr jet in t-channel production at leading and 
next-to-leading order for different scales. 

What stabilizes the pr-ordered jet at NLO if the "d" jet is so unstable? The answer is 
that when the d quark becomes softer, the parton that recoils against the top quark is a 
u that came from initial-state splitting of a gluon. Hence, the kinematically stable jet has 
a different flavor content than the LO jet. We are fortunate that we cannot distinguish 
experimentally the difference between a jet that started from an up or down quark. If 
the jet is produced with a charm or bottom quark, however, we have to track it explicitly 
because it might be taggable. In the case of t-channel single-top-quark production, this is 
CKM-suppressed and occurs less than 1% of the time in the light-quark initiated side of the 
diagram. Once we look at p^-ordered jets, we will have to keep track of how often the leading 
jet includes a b quark from the other incoming hadron. This will not cause any stability 
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problems, however, since there are no overlapping singularities between the light-quark and 
heavy-quark initiated parts of the diagrams. 



A. t-channel NLO distributions 

We have seen the stability of the distribution of the highest-pT jet between leading and 
next-to- leading order. There are two historically relevant questions. First, how well are 
the leading jet and top quark distributions modeled by qb tq' [Fig. Ufa)] or qg tbq' 
[Fig. m^b)] parton-level calculations? Second, how do typical cuts effect the distributions of 
the top quark and jets? The analysis of leading-order qb — > tq' production plus showering 
is postponed until Sec. IIV| when we observe the failure of the current event generators to 
reproduce the jet spectrum. 




FIG. 4: Representative leading-order Feynman diagrams for (a) t-channel production {qh tq'), 
and (b) PF-gluon fusion {qg tbq'). 

The analytic form of the t-channel cross section has the same form as double deep- 
inelastic-scattering (DDIS) jl^. llS^. Color conservation enforces a natural factorization of 
the scales. The fermion line in Fig. EJ^a) that does not include a top quark probes a proton 
with the DIS scale Q^, which is identical to the virtuality of the W boson through NLO. 
The fermion line that does connect to a top quark sees the DIS scale for massive quarks 
of + m^. Another reasonable choice for both scales is the mass of the top quark. The 
difference between using the DDIS scales, or the top-quark mass is less than 3% at next- 
to-leading order. But the cross section and distributions do not change from LO to NLO 
when the DDIS scales are used. This is an indication of the relation of this calculation to 
the extraction of parton distribution functions, and DIS in general. 

Even though the leading-order diagram for t-channel single-top-quark production was 
known to include a b quark in the initial state, the process was modeled by using the tbj 



diagram shown in Fig. Efb) at leading order in Refs. 
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62| . To the extent that the analytic 



form of the cross section is expected to look like double deep-inelastic-scattering [12,1131) this 
diagram was expected to give a reasonable approximation of the shape of the NLO top-quark 
distributions. Indeed, the pxt and pseudorapidity rjt distributions of the top quark are well 
modeled by the LO tbj calculation. These distributions for t production at the Tevatron (a 
1.96 TeV pp collider) are shown in Figs. ElandlHl 
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FIG. 5: Transverse momentum of the top quark pTt at LO, NLO, and from the diagram tbj, all 
normalized to the t-channel NLO cross section. 



The diagram tbj should have no visible effect on the distributions of the highest-pr jet. In 
Figs.[7|and|Hlwe see that there is no difference between the true LO distributions, and the one 
obtained from tbj alone. This is an additional confirmation that there is a real separation 
of physics between the Wtq vertex, and the vertex with no top quark. While this is a 
nice confirmation of the double-DIS interpretation, this means that the distributions for the 
highest-pT jet are more central, and slightly harder than indicated by the old approximation. 
Hence analyses, like that in Ref. slightly underestimated the true signal. 

Finally, let us examine the jet distributions as a first step in understanding the necessity of 
the next-to-leading-order calculations. In particular, we focus on the distributions that will 
be required for Sec. IIVI The simple cuts in Table HI separate the signal into two subsamples: 



t-|-l jet exclusive (Wjj in the experiments, or the "jet veto" of Ref. |2^), and t-|-2 jet inclusive 
{Wjjj in the experiments). Jets are clustered with a /c-r algorithm using AR/^^ < 1.0 (similar 
to a fixed cone of 0.7) and are observed if they have pxj > 15 GeV and \rij\ < 2.5. No cuts 
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FIG. 6: Pseudorapidity of the top quark tjt at LO, NLO, and from the diagram tbj, all normalized 
to the t-channel NLO cross section. 
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FIG. 7: Transverse momentum of the highest-p^ jet pTji at LO, NLO, and from the diagram tbj, 
all normalized to the t-channel NLO cross section. 



are placed on the top quark since its decay is not modeled. 

In Figs. El and ^1 we see the transverse momentum and pseudorapidity distributions for 
the NLO jet in the t + 1 jet final state, distinguished by whether there is a 6 inside the jet. 
The distinction must be made, because the jet containing a 6 is both potentially taggable, 
and has a significantly different shape in the distributions. Notice that non-6 jets exhibit a 
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FIG. 8: Pseudorapidity of the highest-pT jet r]j^ at LO, NLO, and from the diagram thj, all 
normalized to the t-channel NLO cross section. 

TABLE L Cuts on the py-ordered jets that define each signal region for comparison at the Tevatron. 
No cuts are placed on the top quark since it is not decayed. 

|Ai2fe^| < LO |Ai2cone| < 0.74) 
PTj > 15 GeV, \T]j\< 2.5 
Wjj = t + 1 jet ("jet veto") 
Wjjj = t + 2 jets inclusive 

double-peak structure, sometimes referred to as a "forward" jet. Whereas, the additional b 
jets tend to be very central. 

In Figs. HH 1121 CSl andHHwe see the distributions for jets in the t + 2 jet final state. The 
distinction between t plus one and two jet samples is strongly dependent on the cuts. Just 
as in the case of the one jet sample, the b jets are significantly more central. However, they 
actually peak in the opposite hemisphere from the non-6 jets. Approximating the additional 
radiation with just a b jet, as in tbj, greatly exaggerates the centrality of the events. Hence, 
the only way to correctly model the signal is to include all production diagrams in a complete 
next-to-leading-order calculation. 

In theoretical studies that used the tbj diagram to approximate the final state, all ad- 
ditional radiation was implicitly assumed to come from an extra b jet. However, the NLO 
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FIG. 9: Transverse momentum of the highest-py jet p^j^ in t-channel production when it contains 
a b, normahzed to the rate for when it is a non-6 jet in the tj final state. 
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FIG. 10: Pseudorapidity of the highest-py jet rjj-^ in t-channel production when it contains a b, 
normahzed to the rate for when it is a non-6 jet in the tj final state. 

calculation tells us that only about 2/3 of the additional radiation includes a b jet. Fig. 
demonstrates that the contamination from additional b jets in the tj signal is highly depen- 
dent on the choice of cuts. For the loose cuts shown, about 15% of the events have an extra 
b jet, but tighter cuts on the pseudorapidity will greatly enhance the relative contribution 
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FIG. 11: Transverse momentum of the highest-p-r jet PTjx t-channel production when it contains 
a 6, normahzed to the rate for when it is a non-6 jet in the tjj final state. 
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FIG. 12: Pseudorapidity of the highest-py jet r]j-^ in t-channel production when it contains a 6, 
normalized to the rate for when it is a non-6 jet in the tjj final state. 

from events with an extra h jet. 

The next-to-leading-order calculation is necessary to properly account for the absolute 
cross section in each of the four final states: tj, t6, tjj, and thj. In Table HTl we see the 
contribution to each final state for the simple cuts given in Table HI The h jet is the highest- 
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FIG. 13: Transverse momentum of the second highest-pT jet pTj2 in t-channel production when it 
contains a 6, normahzed to the rate for when it is a non-6 jet in the tjj final state. 
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FIG. 14: Pseudorapidity of the second highest-pr jet r]^^ in t-channel production when it contains 
a h, normahzed to the rate for when it is a non-6 jet in the tjj final state. 

Pt jet in 1/3 of the events in the thj final state. The relative cross sections are completely 
insensitive to modest cuts on the top quark itself. However, it bears repeating that the 
percentages will change drastically if the cuts on the jets are changed. Hence, a careful 
matching of the NLO calculation for each choice of cuts is essential to the proper modehng 
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of these events. 



TABLE II: t-channel and s-channel cross sections (pb) for t (or t) production at a 1.96 TeV pp 
collider, and relative percentages of b jets vs. non-6 jets j in the t + 1 jet exclusive and t + 2 jet 
inclusive final states. These numbers are only valid for the cuts given in Tabled Tighter cuts on 
pseudorapidity typically increase the percentage of b jets in the t-channel final states. 

t-channel tj 0.475 pb (87%) tjj 0.076 pb (35%) 

tb 0.073 pb (13%) tbj 0.139 pb (65%) 

s-channel tb 0.310 pb (97%) tbj 0.084 pb 

tj 0.011 pb (3%) 



B. s-channel NLO distributions 



Production of a single top quark 
Drell-Yan production 



igle top quark 



through an s-channel W boson looks very much like 



63|. In Figs. [ISl Cni El and fTHl we see the transverse 
momentum pt and pseudorapidity r] of the top quark and 6 jet. The large mass of the top 
quark protects the distributions from changing shape at next-to-leading order. Hence, it is 
an excellent approximation to use leading-order distributions and a i^'-factor of 1.54 as long 
as we restrict our attention to the tb {Wbb) final state. 

Contributions to the tj and tbj final states arise only at next-to-leading order, and will 
contribute a small correction to the total cross section. In Table ITTl we see that the tj channel 
only contributes 3% to the t + 1 jet final state, and that the ratio between the t + 1 jet and 
t + 2 jet final states is about 4 to 1. The separation between the final states is discussed 
in detail in Sec. IIVBI where we see how the event generators fare at reproducing the jet 
distributions. 



C. EfTect of jet definitions 

Because next-to-leading-order calculations are jet calculations, and not parton calcula- 
tions, results are completely dependent on the choice of jet definition. This means that no 
prediction is meaningful outside of the context of the particular jet definition used. Histor- 
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FIG. 15: Transverse momentum of the top quark pTt in s-channel production at NLO, and LO 
times a ii'-factor of 1.54. 
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FIG. 16: Pseudorapidity of the top quark rjt in s-channel production at NLO, and LO times a 
iT-factor of 1.54. 

ically, jets have been discussed using fixed cones of size Ai?cone < 0.7, where A_Rcone is the 
distance between particles in the pseudorapidity r^-azimuthal angle plane. This definition 
suffers from an infrared instability, because the number of jets can change as particles fluc- 
tuate across the fixed boundary. One infrared safe definition is to choose an algorithm 
that adaptively combines subjets as a function of their relative transverse momenta /ct- For 
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FIG. 17: Transverse momentum of the 6-jet prpi in s-channel production at NLO, and LO times a 
if-f actor of 1.54. 
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FIG. 18: Pseudorapidity of the 6-jet ry^ in s-channel production at NLO, and LO times a K-factor 
of 1.54. 

this study the kr cluster algorithm is used with a merging parameter /S.Rkj, < 1.0, because 
it produces jets that arc very similar to those found when using a fixed-cone of size 0.7. 

In principle, the same jet definition must be used for both the theory and the experimental 
analyses. In practice, the NLO calculation will be matched to an event generator that is 
tuned to the experimental data. The NLO distributions are valid only under the jet definition 
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used to create them. If the experiment decides to use smaller cones, e.g., a cone of size 
AR < 0.4 is currently popular for top physics, the NLO calculation must be reevaluated. 
In Fig. Uni we see the ratio of the distributions pTji{A.R) over ^^^^(0.7), where ji is the 
highest-pT jet in t-channel production, and approximate fixed-cone units label the curves 
(the actual ARk^ used are larger by a factor of 1.35). We see there is a 10% change in 
the slope of the distribution between different typical cone sizes (0.4, 1.0) and the standard 
cone size 0.7. A mismatch between the cone size in the theory calculation and experimental 
calculation could lead to a larger error than all theoretical uncertainties combined. 
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FIG. 19: Ratio da{R)/dpTji to da{R = 0.7) /dpTji for the highest-pr jet in t-channel production. 

When cuts are made on the leading jet, they are typically made between 15-25 GeV. 
From Fig. [7| we see that this is close to the peak of the transverse momentum distribution. 
The jets that are measured have a falling pr spectrum. Hence, Fig. UHl suggests that it is 
advisable to use the largest jet size possible. This would reduce jet-energy corrections for 
out-of-cone effects, and increase the signal acceptance, while having little effect on the Wjj 
backgrounds. The ti background will increase slightly, but is not expected to overcome the 
increased signal. 
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D. Effect of bottom-quark mass 

When there is an additional bottom quark in the t-channel cross section, we might be 
concerned whether its mass is relevant to the calculation. The same almost-coUinear singu- 
larity that must be resummed in the inclusive cross section appears in the distribution of 
this additional jet at low transverse momentum. In Fig. 1201 we see the transverse energy of 
the extra bottom jet Ej'i with no cuts, compared to distribution from leading-order tbj with 
a finite b mass. The curves are very similar (< 10%) above 10 GeV and differ by only a few 
percent by 15 GeV. This is confirms the standard belief that a massive particle behaves the 
same as a massless particle in the large momentum limit. 
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FIG. 20: Effect of the 6-quark mass on the 6-jet transverse energy Erp^ distribution for tbj produc- 
tion with no cuts. 

A more realistic comparison is to consider the difference in jet distributions after cuts. 
The 6 jet with and without the bottom-quark mass included is presented in Figs. |^ and 122 
where a minimal set of cuts for potentially observable jets is defined to have ptj > 5 GeV, 
and \rij\ < 4. In this more realistic case, there is very little difference between including the 
bottom quark mass or not. In any situation where the 6-jets are expected to be taggable, 
the 6-quark mass will have less than a 1% effect on the measured distributions. 

In the s-channel cross section, the extra bottom jet is produced in the final state, and 
is protected from any singularities already at LO by its recoil against the heavy top quark. 
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FIG. 21: Effect of the 6-quark mass on the 6-jet transverse energy Erp^^ distribution for thj produc- 
tion with cuts on both jets pxj > 5 GeV and \7]j\ < 4, and no cut on the top-quark decay products. 
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FIG. 22: Effect of the 6-quark mass on the 5-jet pseudorapidity 7]i distribution for tbj production 
with cuts on both jets pxj > 5 GeV and \rjj\ < 4, and no cut on the top-quark decay products. 

We see in Figs. |2S1 and |^ that the massless and massive b distributions are already indis- 
tinguishable without any cuts. Since there is no visible effect in the measurable portion of 
the distributions, we have ignored the bottom-quark mass throughout. If, in the future, it 
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becomes desirable to measure the third jet below 10 GeV, then mass effects should be taken 
into account. 
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FIG. 23: Effect of the 6-quark mass on the b-jet transverse momentum distribution for s-channel 
single-top-quark production with no cuts. 
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III. UNCERTAINTIES 



Because the measurement of | VJ^p is directly proportional to the measurement of the cross 
section, we must understand the uncertainties in the extracted cross section. Historically, 
theoretical uncertainties have been expected to dominate the cross section measurement. In 
particular there are two significant sources of uncertainty, the uncertainty in the quark-gluon 
luminosity hidden inside the b parton distribution function (PDF), and the modeling of the 
distributions of a third jet that would cause the event to be vetoed. 

In the past, the PDF uncertainty for t(s)-channel single-top-quark production was ap- 
proximated by assuming that it was the same as the gluon(quark) uncertainty. This crude 



approximation was shown to underestimate the uncertainty in the t-channel case 
of 1.5 65|. This study uses the "Modified Tolerance Method" introduced in Refs. 



)y a factor 
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and CTEQ6M PDFs f67|] to provide the first complete estimation for both s-channel and 
t-channel processes. Results appear in Table IIIII The uncertainty at the LHC is much 
smaller than expected. This is partially attributable to the fact that we are probing the 
best measured region of the gluon PDFs, which in turn lowers the b uncertainty. There is 
also an anti-correlation in the quark-6 luminosity that stabilizes the calculation. A more 
conservative estimate would use just the b PDF uncertainty of ^1,2%- Newer CTEQ6M1 
PDFs j3| harden the gluon and sea quarks at moderate fraction of proton momentum x. 
This gives a t-channel cross section for t production of 1.015 pb (a 2.5% increase) at the 
Tevatron, which is well within the current uncertainty estimate. In all cases, there is no 
change in the shapes of the kinematic distributions. 

The shapes of all distributions at the Tevatron and LHC are completely insensitive to 
changes of top-quark mass within a few sigma of its currently measured value. Nevertheless, 
the recent average top-quark mass from run I data of = 178.0 ± 4.3 GeV [GsJ reduces 
the overall normalization of the s-channel cross section for t (or t) production used in this 
paper by 7% to 0.41 pb at a 1.96 TeV pp collider, and the t-channel cross section by 4.8% 
to 0.94 pb. At the LHC (a 14 TeV pp collider), the s-channel cross sections at 178 GeV are 
6.17 pb and 3.85 pb, for t and t production, respectively. The t-channel cross sections are 
152.5 pb and 88.7 pb at the LHC, for t and i production, respectively. Data from run II of 
the Tevatron should reduce the top-quark mass uncertainty to ±2-3 GeV, but it will remain 
the dominant theoretical uncertainty for s-channel production until ~ 30 fb~^ of data are 
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accumulated. 

The complete list of the theoretical uncertainties in the total cross sections are presented 
in Table UTTl for the Tevatron and LHC. The uncertainty due to choice of scale is extensively 
discussed in Ref. and has no effect on the shapes of the distributions. The mass of 
the b quark, and the error in = 0.118 ± 0.003, play an insignificant role in the overall 
uncertainty for measurable quantities. For convenience, the NLO cross sections and their 
current total uncertainties are listed in Table HVl using CTEQ5M1 PDFs [t^. t and t cross 
sections are the same at the Tevatron, but are listed separately at the LHC since it is a pp 
collider. 

TABLE III: Uncertainties of the next-to-leading-order cross section at the Tevatron (pp) and LHC 
(pp). 



Process 


x5mt{GeV) fj,/2-2fi 


PDF 


b mass 




s-channel pp 


-2.330/ 
+2.71/0 


+5.70/ 
-5.0/0 


+4.7q/ 
-3.9/0 


< 0.5% 


±1.4% 


pp 


-1.97 0/ 
+2.26/0 


±2% 


+3.30/ 
-3.9 /O 


< 0.4% 


±1.2% 


t-channel pp 


-1.6 0/ 
+1.75/0 


±4% 


+11.30/ 
-8.1 /O 


< 1% 


±0.01% 


pp 


-0.73q/ 
+0.78/0 


±3% 


+1.3q/ 
-2.2/0 


< 1% 


±0.1% 



Once we begin to calculate exclusive final states, there are additional uncertainties that 
arise from modeling the kinematic distributions. In Sees. Ill Al and IlIBI we saw that NLO 
distributions of the top quark and leading jet were similar to LO distributions including 
i^'-factors if physically motivated scales were chosen. The main uncertainty arises in how 
well the separation between t ± 1 jet and t ± 2 jet distributions is modeled. Since the NLO 
shapes of the jets are stable, the remaining uncertainty is entirely due to the cuts on the 
t ± 2 jet final states, which are only evaluated at LO. 

Typical uncertainties may be derived by observing the change in the normalization of 
the jet distributions when the scales are varied by a factor of 2. The uncertainties for each 
final state depend on the exact cuts that are used. In Table 1 list the uncertainties for 
the cuts in Table E] and cross sections in Table IHI I also provide the complete percentage 
uncertainty in parentheses, assuming 6mt = 4.3 GeV, and the uncertainties in Table IIIII 
The uncertainty for a final state containing a top quark plus one or two jets is presented 
first. The subsamples containing zero or one taggable b jets are also listed. 

23 



TABLE IV: Next-to-leading-order inclusive cross sections for single-top-quark production at the 
Tevatron and LHC for rrit = 175 GeV. Cross sections are evaluated with CTEQ5M1 PDFs, and M^g 
or the DDIS scales {fii = Q^, fJ-h = + ^t)^ s-channel or t-channel, respectively. Theoretical 
uncertainties are added in quadrature, and are based on an uncertainty in the top-quark mass of 
4.3 GeV. 





Vs 

V 


(T1S.TT r\ ( T)Y)] 


s-channel 


1.8 TeV pp (t) 
1.96 TeV pp (t) 
14 TeV pp (t) 
14 TeV pp (i) 


n 077 +0.052 
' -0.045 

n 442 +0 061 

0-00 -0.63 

4 no +0-43 

^■^^ -0.39 


t-channel 


1.8 TeV pp (t) 
1.96 TeV pp (t) 
14 TeV pp (t) 
14 TeV pp {t) 


0.73 t'ofs 
0.99 tin 
155.9 

90.7 tli 



For s-channel production, the tb final state is well modeled, and the uncertainty is com- 
pletely dominated by the top-quark mass measurement. The kinematic uncertainty is less 
than 2.5%. The tj and tbj final states are LO cross sections, and so are more sensitive to 
fluctuations of the additional radiation. The tj and tbj cross sections have approximately 
l^g% kinematic uncertainties in their normalizations, but have very low rates, and are not 
necessary for observation of the signal. 
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In the t-channel, the uncertainty in the tj final state was previously estimated 
±15%, based solely on how often the tbj Feynman diagram produced a b jet with p^b < 15 
GeV. The same large logarithms that lead to the need for the introduction of a 6 PDF, 
also produce a large uncertainty in the shape and normalization of the b jet distributions. 
Changes in the two scales in the problem (Q^ and + mf) tend to accidentally cancel. 
Therefore, an accurate estimation of uncertainty requires varying the scales independently. 

The full NLO calculation tells us that the kinematic uncertainty of the jet in the tj final 
state is only about 6.5% if flavor is ignored. The decrease comes from two effects. First, only 
2/3 of the additional radiation that is cut on includes the poorly modeled b jet distribution. 
Second, the addition of a cut in pseudorapidity removes the forward b jets that are missed. 
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TABLE V: Uncertainties in normalization of jet distributions (pb) in the t + n jet final states for 
t (or t) production at a 1.96 TeV pp collider. These uncertainties are valid for the cuts in Tabled 
and cross sections in Table ITTI Total percentage uncertainties are listed in parentheses, and include 
the uncertainties from Table HTll added in quadrature assuming 5mt = 4.3 GeV. 



Process j flavors t + 1 jet t + 2 jets 

^-channel ah pb pb {t\l%) 

6 jets pb pb (111^ 

16 jet pb IHg^ pb (1^^% 

.-channel ah t'oZ Pb t'of.l pb 

1 I, ;„j. +0.005 f+13o/\ +0.013 „u /+24(v 
11 /O; 



16 jet tHS^pb(l^^ ^ 
6 jets ;Hg^pb(±^^%) 



Hence, the kinematic uncertainty in the tj final state is smaller than either the PDF or 
top-quark mass uncertainties. The NLO distribution of tj when the jet does not contain a 
b is known to 1i[q%- This is still better than, but comparable to the PDF and top-quark 
mass uncertainties. If tighter cuts are made, this uncertainty will improve. Finally, the LO 
tb final state still has an uncertainty of about lti2%- general, more inclusive cross sections 
are better estimated. 

There are two remaining uncertainties that can appear in the theoretical calculation, but 
which are not estimated here. In Sec. II VI a matching between the NLO calculation and event 
generators is proposed. Uncertainties arising from the normalization of the distributions 
produced by the event generators can only be calculated when the fits are performed. The 
last uncertainty comes from the fact that the top quark is not decayed in this analysis. A 
leading-order estimate already limits this uncertainty to less than 2% for the distributions 
and search strategy discussed here. However, if angular correlations are measured 0, , 
or extreme cuts are used, then this uncertainty may be enhanced. A simple remedy would be 

ete published next-to-leading- 



to include the decay of the top quark by combining the comp 
order spin-density matrices for production j3] and decay 



72| into a single program. 
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IV. COMPARISON TO EVENT GENERATORS 



The motivation for performing differential calculations at next-to-leading order is to im- 
prove the connection between theory and experimental results. Typically, this connection 
is made through the use of a showering event generator. These event generators begin with 
leading-order matrix elements, and model the effects of soft initial and final-state showering, 
and hadronization. The essential question is whether these modified leading-order calcula- 
tions are sufficiently predictive for the final states that will be measured. For PYTHIA ^6]]| 
and HERWIG jo^, the answer is a resounding no. 



A. t-channel production in HERWIG and PYTHIA 

We begin by examining t-channel single-top-quark production at the Tevatron (see 
Ref. for a summary of the situation at the LHC). PYTHIA and HERWIG both pro- 
duce events by evaluating the hard cross section represented by Fig. 123 Beginning with 
the correct LO diagram for t + 1 jet production, the programs apply soft showering to the 
final-state quarks (which does not change the primary direction of the overall jet), and use 
strong-angular-ordered showering for incoming quarks. Using the cuts of Table HI we see in 
Figs. Iini and 1221 that the extra jet is reasonably well represented if it is not a b jet. PYTHIA 
and HERWIG have been normalized to the inclusive NLO cross section (the i^-factor is the 
same if you normalize to the t + 1 jet inclusive cross section after cuts). The main discrep- 
ancy here is that PYTHIA overestimates the signal by 8%, while HERWIG underestimates 
it by 5%. 

Unfortunately, both PYTHIA and HERWIG predict completely incorrect distributions 
for the 6 jet in the t+6-jet sample. In Fig.EHlwe see that the transverse momentum predicted 
by the event generators is too soft. While we can imagine trying to harden the spectrum by 
tuning parameters of the event generator, we cannot tune the pseudorapidity distributions 
of Fig. |2ni The showering event generators are producing b jets that are too far forward. 
Once central cuts are placed, PYTHIA (HERWIG) underestimates the number of events in 
this channel by at least a factor of 2.75 (1.4); tighter cuts produce a greater underestimate. 
The origin of the problem is that the showering event generators only produce additional b 
jets from the soft angular-ordered initial-state radiation. All jets produced this way lead to 
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FIG. 25: A symbolic representation of t-channel single-top-quark production as generated by HER- 
WIG and PYTHIA. Soft angular-ordered showering is applied to the hard matrix element for 
ub — ^ td. The radiation is hadronized, and resummed using a jet definition. 
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FIG. 26: Transverse momentum PTqi of the observed non-6 jet in the t + 1 jet final state ("jet 
veto") of t-channel production from NLO (solid), PYTHIA (dashed), and HERWIG (short-dashed) 
normalized to NLO. 



distributions that are too soft and too far forward. 

The failure of PYTHIA and HERWIG to model the t + 6-jet final state has several 
important consequences: Analyses and acceptance studies using PYTHIA or HERWIG will 
lead to approximately a factor of 3 underestimate of the cross section in the Wbb bin. This 

means the cross section for discovery is larger than currently estimated. However, this also 
means that the background to other physics analyses has been severely underestimated. 
The obvious case is for s-channel single-top-quark production, which is distinguished from 
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FIG. 27: Pseudorapidity r/g^ of the observed non-6 jet in the t + 1 jet final state ("jet veto") of t- 
channel production from NLO (solid), PYTHIA (dashed), and HERWIG (short-dashed) normalized 
to NLO. 

t-channel production by the number of h tags. Also, the Higgs analysis {WH Whh) 
will have a significantly larger contribution from t-channel production than estimated. For 
lighter Higgs masses (< 130 GeV) single-top-quark production may be a more significant 
background than ti production. 

The shapes of the distributions in the t + 2-jet final state are also not well-modeled by 
HERWIG and PYTHIA. In Fig. OH we see the transverse momentum and pseudorapidities 
of h and non-6 jets scaled to the NLO event rate. The non-6 leading jets qi are reasonable 
in shape, but either over or underestimated in rate by up to 15%. Other distributions are 
off in rate by up to a factor of 2, and the pseudorapidities are poorly reproduced by both 
HERWIG and PYTHIA. A method for correcting these samples is discussed in Sec. IIV Gl 

B. s-channel production in HERWIG and PYTHIA 

The s-channel cross section is not well represented by either HERWIG or PYTHIA. In 
Figs. l!TTti!?Kl we see that the event generators overestimate the number of events in all final 
states that include additional radiation, and underestimate the number in the t + 6-jet 
iWhh) final state. The origin of this artifact in PYTHIA is that s-channel production is 
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FIG. 28: Transverse momentum pj^;,^ of the observed 6 jet in the t+1 jet final state ("jet veto") of t- 
channel production from NLO (sohd), PYTHIA (dashed), and HERWIG (short-dashed) normaUzed 
to NLO. Inset is on a logarithmic scale. 

implemented by producing a W that is set to have the same mass as a real W boson. The 
problem is that the scale is hard-coded to be the W mass. Hence, this extra radiation 
comes from using too small a scale. The scale should be something on the order of the 
invariant mass of the th pair or at least the top-quark mass. HERWIG also uses a 
scale that is too small, but it uses different scales in the hard cross section for different 
events. For example, gluon-initiated events use an average hard scale around 100 GeV, and 
sometimes use a hard scale less than 1 GeV. This artificially overweights those events which 
have additional radiation. 

The shapes of the transverse momentum spectra are reasonably well modeled by both 
event generators, but additional radiation is too sharply peaked at small pseudorapidity. 
Hence, analyses with tighter cuts than shown in Figs. IHTHH^l will demonstrate a larger dis- 
crepancy. The net effect is that s-channel production in the th-]ei iWhh) final state has 
been underestimated by the event generators by at least a factor of 1.32-1.45; whereas in 
the t + 1 non-6-jet final state, the signal is overestimated by at least a factor of 2-2.4. In 
the t + 2 jet final state, both generators overestimate the number of events by at least a 
factor of 1.4-2.3 depending on the exact final state. Since the discovery channel is Whh, it 
should be easier to find s-channel production than previously estimated by the experiments. 
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FIG. 29: Pseudorapidity rj^-^ of the observed 6 jet 61 in the f + 1 jet final state ("jet veto") of t- 
channel production from NLO (solid), PYTHIA (dashed), and HERWIG (short-dashed) normaUzed 
to NLO. 

However, this means that the s-channel contribution to Wbb as a background to low-mass 
Higgs production has been underestimated by the same amount. 



C. Correcting the event generators 

The first general-purpose next-to-leading-order showering event generators for hadron 
colliders are beginning to appear 0|. However, they continue to rely on the ability to 
cleanly separate initial-state and final-state radiation. Single-top-quark production through 
the s-channel exhibits such a separation, but has not yet been added to the programs. 
A corrected leading-order generator using a K-factor is expected to perform well enough. 
On the other hand, there is no clear separation in t-channel production between initial-state 
and final-state radiation. Given the similarity to deep-inelastic scattering, it may be possible 
to construct a next-to- leading-order event generator for t-channel production jr^], but no 
working implementation exists. Nevertheless, in order to compute detector efficiencies it 
is necessary to have predictions at the level of detectable particles (pions, photons, etc.). 
Therefore, I propose a simple matching prescription that bypasses many of the problems 

0. 
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FIG. 30: Transverse momentum px and pseudorapidity r/ for the p-r-ordered jets ji in the t + 2-jet 
final state of t-channel production from NLO (sohd), PYTHIA (dashed), and HERWIG (short- 
dashed), ji and 22 are divided into non-6 jets qi and b jets 6j. Curves are normahzed to NLO in 
each column. 

The goal is to produce a set of events with the proper shape and normalization. The 
most important thing to match is the shape of the distributions in a given final state. Once 
the shapes are matched, generated events may be multiplied by an overall probability that 
is normalized by the NLO calculations after cuts. This prescription is based on the essential 
understanding that the next-to-leading-order calculations are jet calculations. Therefore, 
matching must be done at the jet level. This is accomplished in four steps. 

1. Generate events with MadEvent [76^ or CompHEP T^l using the tree- level Feynman 
diagrams that produce all measurable final states. The t-channel parton-level processes 
for top-quark production are qb tq', qb tq'g, gb tqq', and qg tq'b, where q 
and q' stand for the five light quarks. The s-channel parton-level processes are qq' tb, 
qq' tbg, qg tbq', and gq' — > tbq. Use the DDIS scales {fii = , fih = + 'm'i) 
for t-channel production, and Mtb for s-channel production to minimize changes in the 
shape between LO and NLO. 

2. Feed these events into HERWIG or PYTHIA with full showering turned on. 

3. Create four samples (tj, tb, tjj, and tbj) for t-channel production with minimal cuts. 



31 



b 
-a 



4.5 
4.0 F 
3.5 
3.0 F- 
2.5 r 
2.0 r 
1.5 r 
1.0 r 
0.5 E- 
0.0 



NLO 

1.32xPYTHIA 
1.45xHERWIG 



7j,J <2.5 




L±J I I I I I I I I I I I I I I I I L 



50 100 150 

PTb, (GeV) 



200 



FIG. 31: Transverse momentum of the observed 6 jet in the t+1 jet final state ("jet veto") of s- 
channel production from NLO (soHd) , P YTHIA (dashed) , and HERWIG (short-dashed) normahzed 
to NLO. 
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FIG. 32: Pseudorapidity rjb^ of the observed b jet bi in the t + 1 jet final state ("jet veto") of s- 
channel production from NLO (solid), PYTHIA (dashed), and HERWIG (short-dashed) normalized 
to NLO. 

and a well-defined jet definition, e.g., kr clustering with AR < 1, applied to stable 
particles in the event record (tt^, 7, e^, //^, etc.). s-channel production will have three 
samples {tb, tj, and tbj). 
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FIG. 33: Transverse momentum PTqi of the observed non-6 jet in the t + 1 jet final state ("jet 
veto") of s-channel production from NLO (soUd), PYTHIA (dashed), and HERWIG (short-dashed) 
normahzed to NLO. 



ri 0.5 Ix PYTHIA ---J 
|J' b.41xHERWIG - 




-3-2-10123 



FIG. 34: Pseudorapidity r/qj of the observed non-6 jet in the i -|- 1 jet final state ("jet veto") of s- 
channel production from NLO (solid), PYTHIA (dashed), and HERWIG (short-dashed) normalized 
to NLO. 

4. Produce the same final states with the NLO calculation, and use the same cuts and jet 
definition. Normalize the weights of the events in each of the event generator samples 
to the cross sections of these NLO samples. Remember that the NLO jets should be 
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FIG. 35: Transverse momentum px and pseudorapidity r/ for the p-r-ordered jets ji in the t + 2-jet 
final state of s-channel production from NLO (solid), PYTHIA (dashed), and HERWIG (short- 
dashed), ji and j2 are divided into h jets hi and non-6 jets qi. Curves are normalized to NLO in 
each column. 

i?T-ordered, but that explicit heavy-flavors in the final state must be tracked. 

In principle we would want to apply the jet definition to the event generators at the point 
of clustering. However, this is not possible in PYTHIA because of the order in which events 
are processed. Furthermore, both event generators do some slight reshuffling and energy 
rebalancing at a later point in the process. The NLO calculation is not aware of long-range 
final-state effects, but constraints in the event generators on the direction of particles coming 
from the clustering stage maintain the correspondence with the perturbative calculation. 
Matching to jets produced out of stable particles has been found to work to better than 
1% (the limit of the statistical accuracy checked) for the simpler case of s-channel W 
production j^^l^^, where the event generators perform well. In general, a systematic study 
of the uncertainty in the matching should be performed for each set of cuts. 



V. CONCLUSIONS 



Precise measurements of single-top-quark production in both the s-channel and t-channel 
modes will begin to appear within the next few years using data from run II of the Fermilab 
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Tevatron, and eventually the CERN LHC These cross sections are distinguished by the 
frequency of b jets occurring in the + 2 jet sample, with Wbb (Wbj) most often coming 
form s-(t-)channel production. These measurements will provide direct extraction of the 
CKM matrix element Vtb, and may allow discovery and interpretation of new physics. 

Historically, the ability to interpret the measurement of the single-top-quark cross sections 
was expected to be dominated by theoretical uncertainties. In Sec. IIIH we see that these 
uncertainties, while still large, are under much better control. The uncertainty in the top- 
quark mass completely dominates the prediction for the s-channel cross section, but this 
will improve with additional data. The measurement of the s-channel cross section will be 
statistics limited at the Tevatron, and difficult to distinguish from backgrounds at the LHC 
I23I 1. Hence, the theoretical uncertainty in s-channel production should be smaller than the 
expected experimental accuracy at either a high- luminosity Tevatron, or the LHC. 

Current theoretical uncertainties in t-channel production limit extraction of Vtb to about 
±8%. The next-to-leading-order distributions presented in Sees. Ill Al and IIV Al improve 
the kinematic uncertainty in t-channel production, so that parton distribution function 
uncertainties will dominate our ability to interpret the data. This uncertainty could be 
reduced by measuring the ratio of t-channel production to Zb production. However, the 
latter cross section is small, and thus we should consider alternate possibilities, such as 
Wc or Wj exclusive final states. Any improvement in the measurement of the b or gluon 
distributions at scales near the top-quark mass, and proton momentum fraction x near 0.1, 
will directly reduce the remaining theoretical uncertainty. 

The availability of fully differential next-to-leading-order calculations for hadronic pro- 
cesses has begun a new era of precision comparisons between theory and experiment. In 
order to fully realize the potential of these theoretical improvements, we must first under- 
stand how to interpret the calculations. In Sec. Ull we see that NLO calculations produce jet 
distributions in the final state, rather than parton distributions. This forces us to describe 
the results in terms of i?T-ordered jets, rather than NLO distributions of a specific quark. It 
further ties the theoretical results to the explicit choice of jet definition. We see in Sec. Ill C\ 
that the effect of jet definition on the shapes of distributions is larger than all theoretical 
uncertainties combined. 

Finally, fixed-order calculations are generally compared to showering event generators, 
such as PYTHIA and HERWIG, rather than directly to data, in order to understand detector 
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efficiencies. Therefore, we must understand how to compare next-to-leading-order jets to 
jets produced from showering. A simple matching scheme is presented in Sec. IIV CI that uses 
the NLO calculation to normalize the measurable event samples. Ultimately, we would like 
to have a next-to-leading-order event generator that can handle initial-state and final-state 
interferences. 

Now that the theoretical predictions for single-top-quark production are under control, 
a few questions remain regarding the decay of the top quark. One of the most interesting 



features of sinj 
top quark ^ 



e-ton 
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quark production is the possibility of measuring the polarization of the 



22] • The program used to produce distributions for this paper, ZTOP, 
does not decay the top quark. However the complete analytic NLO spin-density matrices 
for production do appear in Ref. and could be combined with the known NLO decays 
fromRef. 3- 

A spin-dependent program would allow an estimate of the uncertainties associated with 
using alternate search strategies, such as cutting on Ht, the scalar sum of the transverse 
energies of the particles in the event. It is clear from the distributions in Sec. IIV| that the 
Ht for single-top-quark production is larger than predicted at LO in the -|- 2 jet inclusive 
sample. The inclusion of top-quark decay is necessary to determine whether this is a serious 
problem or not. Ultimately, a program that includes decays will be necessary to determine 
whether the single-top-quark background to new physics, such as Higgs production, has been 
greatly underestimated, or whether specific cuts accidentally reduce the discrepancies. For 
now, however, we have a reliable theoretical calculation, and method for matching to event 
generators, that will allow for accurate interpretation of single-top-quark production data 
when it is discovered. 
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